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AN INTEGRO-PDE MODEL 
FOR EVOLUTION OF RANDOM DISPERSAL 

KING-YEUNG LAM AND YUAN LOU 


Abstract. We consider an integro-PDE model for a population structured by 
the spatial variables and a trait variable which is the diffusion rate. Gompeti- 
tion for resource is local in spatial variables, but nonlocal in the trait variable. 
We focus on the asymptotic profile of positive steady state solutions. Our re¬ 
sult shows that in the limit of small mutation rate, the solution remains regular 
in the spatial variables and yet concentrates in the trait variable and forms a 
Dirac mass supported at the lowest diffusion rate. Hastings and Dockery et 
al. showed that for two competing species in spatially heterogeneous but tem¬ 
porally constant environment, the slower diffuser always prevails, if all other 
things are held equal I13II15I . Our result suggests that their findings may well 
hold for arbitrarily many or even a continuum of traits. 


1. Introduction 

In this paper, we focus on the concentration phenomena in a mutation-selection 
model for the evolution of random dispersal in a bounded, spatially heterogeneous 
and temporally constant environment. This model concerns a population structured 
simultaneously by a spatial variable x € D and the motility trait a € A of the 
species. Here D is a bounded open domain in and A = [a, a], with a > a > 0, 
denotes a bounded set of phenotypic traits. We assume that the spatial diffusion 
rate is parameterized by the variable a, while mutation is modeled by a diffusion 
process with constant rate > 0. Each individual is in competition for resources 
with all other individuals at the same spatial location. Denoting by u(t,x,a) the 
population density of the species with trait a G A at location x G D and time 
t > 0, the model is given as 

{ Ut = aAu + [m{x) — u(x, t))] u -I- e^Uaa, x G D,a G {a, d), t > 0, 

= 0, X G dD, a G (a, a), t > 0, 

Ua = 0, X G D,a G {a, a} A > 0, 

u(0, X, a) = uo{x, a), x G D, a G {a, a). 

A = denotes the Laplace operator in the spatial variables, 

pa 

u{x,t) ■= / u{t,x,a)da, 

J 

n denotes the outward unit normal vector on the boundary dD of the spatial domain 
72, and ^ = n ■ V. The function m{x) represents the quality of the habitat, which 
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is assumed to be non-constant in x to reflect that the environment is spatially 
heterogeneous but temporally constant. 

The model dm can be viewed as a continuum (in trait) version of the follow¬ 
ing mutation-selection model considered by Dockery et al. m, concerning the 
competition of k species with different dispersal rates but otherwise identical: 


( 1 . 2 ) 


dt' 


= -I- 




dn 


U^{x, 0) = Uifi{x) 


Ui + J2'j=i 

in D X (0, oo), i = 1,..., fc, 
on dD X (0,oo),i = 1, 
iT\ D,i = 1,..., k, 


where 0 < oi < 02 < ... < ccfc are constants, m{x) G C^{D) is non-constant, is 
an irreducible real k x k matrix that models the mutation process so that Mu < 0 
for all i, and Mij > 0 for i ^ j and > 0 is the mutation rate. 

Model (O was introduced to address the question of evolution of random dis¬ 
persal. In the case when there is no mutation, i.e. e = 0, this question was 
considered in where it was shown that in a competition model of two species 
with different diffusion rates but otherwise identical, a rare competitor can invade 
the resident species if and only if the rare species is the slower diffuser. Dockery et 
al. [13] generalized the work of Hastings m to k species situation, and proved that 
no two species can coexist at equilibrium, i.e. the set of non-trivial, non-negative 
steady states of the system (1121) is given by 


{( 0 „ 1 , 0 ,..., 0 ), ( 0 , 0 ,..., 0 ),..., ( 0 ,..., 


where da is the unique positive solution of 

aA6 + 0{m — 0) =Q in D, = 0 on dD. 

Moreover, among the non-trivial steady states, only (0aj,...,O), the steady state 
where the slowest diffuser survives, is stable and the rest of the steady states are all 
unstable. Furthermore, when k = 2, the steady state (0 q,j,O) is globally asymptot¬ 
ically stable among all non-negative, non-trivial solutions. Whether such a result 
holds for three or more species remains an interesting and important open question. 

Dockery et al. m further inquired the effect of small mutation. More precisely, 
when 0 < e <C 1, it is shown that (HU has a unique steady state U = {ui,U 2 , 
in the space of non-trivial, non-negative functions, such that Ui > 0 for all i, and 
U —>■ {0ai, 0,..., 0) as e —0; i.e. the system (11.21) equilibrates only when the slowest 
species is dominant and all other species remain at low densities. 

It is natural, then, to inquire if the situation in the discrete (in trait) framework 
carries over to the continuum framework. The aim of this paper is to study the 
asymptotic behavior of steady state(s) of (11.11) . Let be any positive steady state 
of (|l.ip . we will show that, as e —>■ 0, 

Ue{x, a) —>■ 5{a — a)0a{x), 

i.e. Ue converges to a Dirac mass supported at the lowest possible trait value a. 
See Theorem 12.31 for precise descriptions of our main results. 

Mutation-selection models for a continuum of trait values have been studied 
extensively, when the phenotypic trait is associated only with growth advantages 
[HlHlEKnKlTllIllEI]. See also [16] for a pure selection model. The consideration 
of a spatial trait is more recent [DEllzlEo]. 
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System dD is also considered in an unbounded spatial domain x S R. A formal 
argument concerning the existence of an “accelerating wave” is presented in [6], 
which provides a theoretical explanation of the accelerating invasion front of cane 
toads in Australia |23) . Rigorous results are obtained when a £ A = [a, a] more 
recently in [Sill]. It can be summarized that the highest diffusion rate is selected 
when the underlying spatial domain is unbounded, which stands in contrast to the 
case of bounded spatial domains we consider in this paper, where the lowest possible 
diffusion rate is selected. 

The rest of the paper is organized as follows: The main results are stated in 
Section!^ Section [3] concerns various estimates on steady states of (11.11) . In Section 
|3]we introduce an auxiliary eigenvalue problem and a transformed problem of (|2.1I) . 
The limit of tie is determined in Section [SJ In Section [51 we analyze the qualitative 
properties of solutions to the transformed problem. The proof of our main result is 
given in Section 0 Finally, the Appendices lAl to ICl establish the existence results, 
the smooth dependence of principal eigenvalue on coefficients as well as a Liouville- 
type results concerning positive harmonic functions on cylinder domains. 


2. Main Results 


In this paper, we consider the asymptotic behavior of positive steady states of 
dm), denoted by u^. That is, satisfies the following mutation-selection equation 
of a randomly diffusing population: 

aAue + (^{ue)aa + [m{x) — ?}.e(x)] We = 0 in Q := D X (a, a), 


( 2 . 1 ) 


dUe 


= 0 


dn 

('^e)a — 0 


on dD X {a, a), 
in D X {a, a}, 


where 

( 2 . 2 ) 


pOc 

u^{x) = / Ue{x,a)da. 
Throughout this paper, we assume 
(A) m{x) is a non-constant function in C{D) such that 


m{x) dx > 0. 


ID 


The existence of positive solutions to dm can be stated as follows: 

Theorem 2.1. Suppose (A) holds, then (12.11) has at least one positive solution for 
all e > 0. 

We postpone the proof of Theorem l2.1l to Appendix lAl For the rest of the paper 
we will focus on the asymptotic behavior of positive solutions of (j2.ip as e —>■ 0. To 
this end, we define the following quantities: 


Definition 2.2. 

(2.3) 


(i) Let 0 q(x) be the unique positive solution of 

aA0 + 0{m{x) — 0) = Q in D, 

= 0 on dD. 

on 


(2.4) 


(ii) For each a € [a,a], we denote the principal eigenvalue and principal posi¬ 
tive eigenfunction of the following problem by cr*(a) and ip*{x,a), respec¬ 
tively: 

aAtjj + (m(x) — 0a{x))'ili + crip = 0 in D, 

= 0 on dD, and -0^ dx = 0^ dx. 
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(Note that by (i), 9ai.x) is a positive eigenfunction for (12.41) when a = a. By 
uniqueness of the (normalized) principal eigenfunction, we have cr*(a) = 0 , 
and ip*{x,(£) = da{x) for x G D.) 

(hi) Denote by 77 * (s) the unique positive solution to 


(2.5) 


J 77 " + (ao — 015)77 = 0 for s > 0 , 

( 77'(0) = 0 = 77 (-|-oo) and 77 ( 5 ) ds = 1, 

where 00,01 are positive constants determined by oi = ^^(o) and oq = 
(ai)^/^Ao, where Aq is the absolute value of the first negative zero of the 
derivative of the Airy function. 


When m(x) = 1, one can easily show that Ue = l/(a —a), i.e. there is no 
selection in the trait variable. Our main result shows that the outcome changes 
drastically when m(x) is non-constant. In fact, Uc concentrates at the lowest value 
in the trait variable, as e —>■ 0. This phenomenon is also known as spatial sorting. 


Theorem 2.3. Let Ue be any positive solution of m- Then for all /3 > 0, there 
exists C > 0 independent of e > 0 such that 

(2.6) Ue{x, a) < exp /3(a — 

in fl = D X (a, a). Moreover, as e —>■ 0 

(2.7) e^/^Ue{x, a) - 9a{x)r]* 

where 9a{x) and rf(s) are given as above. In particular, we have 

pa 

( 2 . 8 ) Ue{x) = / Uf^{x,a) da ^ 9a{x) as e —>■ 0. 

J 


L°°(n) 


As the proof of Theorem 12.31 is fairly technical, we briefly outline the main 
ingredients for readers. Our idea is to establish the “separation of variables” formula 
(EJl) for Tie, by introducing the scaling s = (a — a)/e^/^ and writing Ue{x,a) = 
ipe{x,a)we(x, s), where 7/>e is the principal eigenfunction of —aAip + (ue — m)%f = 
(T'0, subject to the zero Neumann boundary condition and the integral constraint 
= Jjj9^. The main body of our paper is devoted to the proof of following two 
things: (i) As e —>■ 0, Ug — 9g_ uniformly, which implies that 'ife{x,a + e^/^s) —>■ 9a. 
The concentration phenomenon of Ug on the subset D x {a} of D, i.e. (12.81) . is 
established in Section 5 with the help of several key estimates proved in Sections 
3 and 4, such as the L°° estimate of Ug, as well as the limit limg_>.o('ag — m) being 
non-constant; (ii) As e —>■ 0, e~‘^^^We{x, s) converges to 77 ( 3 ) uniformly for some 
function rj which is independent of the spatial variable x. This is done in Section 6 
along with some key estimates established in the earlier sections. 


Remark 2.4. After this work is completed, the authors learned that a closely re¬ 
lated result, under a slightly different formulation, is independently proved by B. 
Perthame and P.E. Souganidis under a different approach, where an intermediate 
trait attains the minimum diffusion rate and an interior Dirac mass is found when 
the mutation rate tends to zero. Apart from the distinction in our approaches, 
we note the following distinct features of our work: (i) A boundary concentration 
is found in our set-up, instead of an interior concentration in |22] which predicts 
different scalings in powers of e; (ii) Our method does not assume the convexity 
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of spatial domain D] (iii) Various detailed L°° estimates and asymptotic limits 
are obtained (Theorem 12.3|) which paves the way to the proof of asymptotic sta¬ 
bility and uniqueness of in a future paper; (iv) The key estimate of the limit 
/io(a;) = lim£_>o ne(a;) — m{x) being non-constant iLemma 13.41) reflects the effect of 
spatial heterogeneity, the underlying mathematical reason for the selection of small 
diffusion rate. See also Proposition 13.71 which makes the connection to [22l Lemma 
4.3]. 


3. Properties of Ue 

In this section we establish various properties of Ue- Recall that is defined in 

( 1221 ). 


Lemma 3.1. There exists some positive constant (5i = Si{a,a,m) independent of 
e such that 

<5i < Ue{x) < 1/^1 in D 

for all e > 0. In particular, 

(3.1) he{x) := Ue{x) — m(x) 
is bounded uniformly in L°°{D). 

Proof. The idea of the upper bound follows from [24]. Define 

(3.2) Ve{x) = / aue{x,a)da. 

J ^ 

Then we have 


(3.3) aue{x) < Ve{x) < aue{x) in D. 

Integrating (EH) over a gives 


(3.4) 


Ave{x) + {m{x) — Ue(x))ue{x) =0 in £>, 

= 0 on dD. 

on 


Let max|) Pe = Ve(xo), then Ue(xo) < ni(xo) < maxjjm (see [HI Proposition 2.2]), 
and by (13.31) . 


a max Me < maxve = Ve(xo) < aue(xo) < amaxm. 

D D D 


Hence we deduce that l[Mel|Loo(_D) and ||ft-e||L~(_D) are bounded uniformly in e, where 
he{x) = Ue{x) — m{x) is given in (13.1|) . 

Next, we show the lower bound of Ue- By (13.31) . we deduce that 


Ue{x) = ke(x)Ve(x) 

for some fce(x) S L°°{D) such that a~^ < ke{x) < a~^. So that Ve is a positive 
solution of 


dv 

—Ave + he{x)ke{x)ve = 0 in D, and ® on dD, 

where we have already shown that h^ = Ue — m is uniformly bounded (in L°°{D)) 
in e. Therefore, the Harnack inequality applies so that 


(3.5) 


max Ve < min 

D D 
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for some constant C > 1 independent of e. Combining with (I3.3p . we have 
(3.6) amaxue < maxiig < C' minve < C^dminfte. 

D D D D 


Now, if we divide (EB by Me and integrate by parts over V, = D x {a, a), we 
obtain 

(3.7) {a-a) ( {u,-m)dx= [ K{x) dadx = [ > q, 

Jd Jn Jn uf 

We deduce by (13.61) and (13.71) that 


C'a ^ ^ 1 /■ .- e N , 1 

- mm Me > max Me > / uJx) dx > t-— r 

a D D \D\ Jd \D\ , 

This establishes the uniform lower bound of Me. 


m{x) dx > 0. 


ID 


□ 


Remark 3.2. Since ||Me||Loo(£)) and ||Me||Ltxj(£)) are bounded uniformly in e, applying 
elliptic LP estimate to (13.41) implies that ||Me||iv2,p(£)) is bounded uniformly in e. In 
particular, there exists sequence —>■ 0 such that converges uniformly on D. 

Lemma 3.3. There exists a constant C > 0 such that for any positive solution Me 
of dHH), 

sup Me < Ce~^. 

D X (a ,a) 

Proof. Choose Xe and such that the supremum of Me is attained at (Xe, a^). Next, 
let C/e(x, r) = Me(x, Oe + cu), then (extending Me to D x [a — eo,a + cq] by reflection 
across the boundary portions D x {a, a} if necessary) one may observe that Ue 
satisfies a uniformly elliptic equation with uniformly bounded (in L°“) coefficients 

f Q;Aa;{7e + t7e,Tr + /le(x)t/e = 0 inL>x[-2,2], 

1^=0 on5Z7x [-2,2 ], 

where a = tte + er is always bounded between [a — cq, a, +eo] C (0, +oo). Hence, 
we may apply the Harnack’s inequality to yield a positive constant C independent 
of e such that 

Ue{x,ae + er) > CUe{Xe,ae) = C sup Me 

Dx (a,a) 

for all X G £>, T G [—1,1]. Hence, 

poi 

Ue{x) = / Ue{x,a)da>Ce sup Me 

Dx{a,a) 

for all e sufficiently small. By Lemma [3.11 we deduce that 

sup Me < C"e“^ 

Dx{a,a) 

for some positive constant C. □ 

By Lemma [3.11 he is bounded in L°°{D) uniformly in e. Therefore, up to sub¬ 
sequences Cj —>■ 0, hej converges weakly in LP{D) for all p > 1. We first prove an 
important property of any subsequential limit h^. 

Lemma 3.4. Let ho be a weak (subsequential) limit of he{x) in LP{D) (p > 1) as 
e —>■ 0, then ho{x) is non-constant in D. 

Proof. Suppose to the contrary that for some c G R, he{x) c weakly in LP{D) 
for all p > 1. 
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Claim 3.5. c = 0. 

By taking e —?> 0 in (1X71) . we deduce that c > 0; i.e. for some c > 0, Ue = 
he + m{x) ^ c + m{x) in U’ for all p > \. 

Next, integrating with respect to a and then x, we obtain 



so that 



where the first inequality follows from (13.81) . and the second inequality from ex¬ 
panding Jjj[ue{x) — {m{x) -|- c)]^ > 0 as 



Hence, by (13.91) . we have c = 0; i.e. he{x) 0 weakly in LP{D) for all p > 1. 

Note that we are done if m < 0 somewhere, since then he{x) > —m(x) > 0 in 
some open subset of D independent of e, which contradicts he 0. For the general 
case of m{x) being possibly non-negative, we continue via a blow-up argument. Let 



It is enough to show that 
Claim 3.6. Ce \ 1 as e —0. 

Assuming Claim [SHI then by definition of Ce, 

Ue{x, a) < CeUe{y, a) for all x,y € D and a < a <a. 

This gives, upon integrating over a € {a, a), 

sup tie(x) < Ce inf Ue{y). 

x€D yeD 

Hence Ue(x) converges to a constant. But this also means that he = Ue — m{x) con¬ 
verges to a non-constant function, as m(x) is non-constant. This is a contradiction. 

It remains to prove Claim 13.61 Assume to the contrary that there exist some 
constant cq > 1, and sequences Cfc —>■ 0, —>■ ao: J/fc G L) such that 

(3.10) Uefc(Xfc,afc) > CoMefc(yfe,afc)- 

Extend Ueto Dx [a —cq, a-l-eo] for some fixed eo small by reflection on the boundary 
D X {a, a}, and define 



Then (13.101) says that for some cq > 1 independent of k 


(3.11) inf Uk{x,0) < — for all k. 



Moreover, Uk satisfies 




). 


Ukix,s)>0 inl^x sup,g^t/fc(x,0) = l. 
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Since afc —>■ ao S [a, a], Ck —>■ 0, the domain of Uk converges to 13 x R as fc —>• 
oo. Moveover, by the uniform boundedness of in k fLemma 13.11) . we 

have for each M > 0 the coefficients of the equation of Uk{x,s) are bounded in 
L°°{D X [—M, M]) uniformly in k. Since sup^g^, Uk{x, 0) = 1, together with p.lll) 
we may apply Harnack inequality to obtain a constant C = C{M) independent of 
k such that 

C~^ < Uk{x, s) < C ioT X G D and |s| < M. 

By LP estimates (applied to D x [—M,M] for each M), there is a subsequence 
Uki that converges uniformly in compact subsets of l) x R. to a positive solution 
of AxUq + {Uo)ss = 0 on 13 X R. (The limiting domain is 13 x R as — —>■ 
—oo and —>■ oo.) Now, we apply Proposition 1C.II for positive harmonic 

functions on a cylinder domain, so that Uq = ci for some positive constant ci. 
Since sup^^]jUk{x,0) = 1 for all fc, we have ci = 1. In particular, we set s = 0 
and find a subsequence Ufc. (a;, 0) converges to 1 uniformly for x G D. This is in 
contradiction to ()3.11|1 and proves Claim 13.61 This completes the proof. □ 

The following result generalizes a key estimate of [22] , proved wherein via Bern¬ 
stein’s method under the additional assumption that D is convex. Although not 
needed for the rest of the paper. Proposition 13.71 enables one to follow the elegant 
Hamilton-Jacobi approach as in [22] to show the concentration phenomenon. 

Proposition 3.7. Let he a positive solution of eu). Then there exists a eon- 
stant C > 0 independent of e such that 




V x'^e 



Ue 

L°°{Q) 

Ue 


Proof. Extend the definition of Ue to 13 x (2a — a, 2a — a) by reflecting along the 
boundary portions D x {a, a}. For each ao G [a, a], define 

Ue{x, t) := Ue{x, ao + er). 

Then is a positive solution to 

A{t, e)AxUe + Ue,TT - he{x)Ue = 0 

in 13 X (e“^(2a — a — ao), e“^(2a — a — ao)) and satisfies the Neumann boundary 
condition on dD x (e“^(2a —a — ao), e“^(2a — a — ao)). Here A(t, e) is a contin¬ 
uous function such that a < A{t) < a. This, together with the boundedness of 
||^e||L“(D) fLemma l3.1L one may apply the Harnack inequality to Z3 x (—1,1) and 
deduce the following. 

Claim 3.8. There exists C > 0 independent of ao G [a, a] and e sueh that 
sup Ue{x,T)<C inf Ue{x,T). 

Dx(-l.l) Dx(-l.l) 

Next, we apply elliptic estimates to C/g in Z3 x (—1,1), so that 
(3.12) 

sup [|t7e,r(a;,0)|-b |V 3 ,C/e(a;, 0 )|] < < C sup Ue{x,T). 

xGD y ’ Dx(-l.l) 

In view of Claim [3^ we deduce for any x G D, 

\UeAx,0)\ + \V^U^ix,0)\ <C inf [/, < CC/,(a;,0). 

Dx(-l,l) 
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i.e. e |Me,c((a;, ao)| + |VxMe(a;, Q!o)| < Cue{x,ao) ioi aWx € D. Since C is independent 
of X, ao and e, this proves Proposition 13.71 □ 


4. Two EIGENVALUE PROBLEMS 


4.1. An Auxiliary Eigenvalue Problem. Consider, for each a > 0 and e > 0 
the eigenvalue problem (recall h^{x) = u^{x) — m{x)) 


(4.1) 


—aAtjj + = aip in 

= 0 on dD and '0^ dx = dx, 


and denote the principal eigenvalue and positive eigenfunction by a^^a) and 0')j 
respectively. At this point, we have not shown how the two eigenvalue problems 
631) and (1^ are related yet. 

For each e > 0, fTe(a) is a smooth function of a > 0 fProposition IB.ll iill. and it 
has a Taylor expansion at a = a: 


(4.2) (Te{a) = (To,e + cri,e(a - o) + cr 2 ,e(a - af + 0{{a - of), 


where cro,e 


tTe(a) and Uk 



o'e(a)- 


Lemma 4.1. Let and be given as above. 

afc 

(i) For each k>0, is bounded uniformly in e > 0 and a G [fi, a]- 

(ii) For eaeh k > 0 and p > 1, •, ci) *5 bounded in W'^’P{D) (and hence 

C{D)) uniformly in e > 0 and a G [a,a]. 

(iii) There exists cq > 0 such that 


lim inf ^ (a) > cn > 0 
e-i-o da 


for all a G [a, a]. 


(iv) 


d(7 

In particular, lim inf tri e = lim inf ^(a) > 0. 
e^O ’ e^O da 

There exist positive constants ci, C 2 such that for all e > 0, 


Cl < 'f’eix, a) < C 2 for all x G D and a G [a, a]. 


Corollary 4.2. There exists C > 0 independent of e such that 


'0£,a 


+ 

L°° {D X (a,c?)) 


'4^e,aa 


L°^(DX (a,Q:)) 


< c. 


Proof of Lemma EH By the uniform boundedness of ||h£||/,cx,(£)) in e (Lemma l3.ll) . 
assertions (i) and (ii) follow from Proposition IB.Sl ib To show (iii), it suffices to 
show, given any sequence Cj —>■ 0, and aj —>■ ao G [a, a], liminfj_>.oo > 0. 

By Lemma lSTl we may assume without loss of generality that for some hg G L°° (D), 
he ho weakly in Lp{D) for all p > 1. Then, in the notation of Appendix iBl 
Proposition IB.5l iii implies that 

d d d 

^ -^Xiiao,ho). 

Since hg is non-constant (Lemma 13.41) . Proposition IB.ll implies that the last expres¬ 
sion is positive. This proves (iii). 

For (iv), suppose that along a sequence Cj —>■ 0 and aj ^ ag > 0, either 
iniioipejix,aj) —^ 0 or sup^, (x, a^) —>■ oo. By the uniform boundedness of 
ll^e||L”(D) (Lemma 13.II) . we may assume without loss that he- converges weakly in 
LP{D) for all p > 1. Hence by Proposition IB.Sl iil. i(ej = Ti{' cij, (ie„ ) converges to 
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Lpi{ ■; ao, ho) uniformly in Z?, and the latter is a strictly positive function in C{D). 
This is a contradiction, and proves (iv). □ 


4.2. A Transformed Problem. By the fact that Ue{x,a) is the principal eigen¬ 
function, with zero as the corresponding principal eigenvalue, of the problem 

—aAcj) — e'^(f>aa + = 0 in = £> x (a, a), 


(4.3) 


■^cj) = 0 on dD X (a, a), and (j)a = Q on D x {a, a}, 


we have the following variational characterization 


(4.4) 

0 = inf Jg [(j)] , 

<#> e H^(n) 

/fi = 1 

where 


(4.5) 

Je[<jl]= [ [a\Vx(l)\‘^ + . 

Jn 

Dehne 


(4.6) 

Se = (a - a)/e^/^. 

and 


(4.7) We{x,s) 

:= Ue{x,a + e^/^s)/'0g(a;, a + e^^^s) for x £ D, 0 < 


where is given by m- Then We satisfies 

-oVa; • {illVxWe) - e'^/^{^'^We,s)s + l/’? Cre{a) - ^ 

in Z? X (0, Se) 


We = 0 


(4.8) 


^We=0 onaHx(0,Se), 


w. 


:.s = -e2/3^We On I? X {0, Se}- 


The corresponding variational characterization can be written as 

M4>] 


(4.9) 

where 

J^[(j)] = 


U^0,e — (tl) — 


inf 


ID Jo 
e 


<l,eH^(Dx( 0 ,s,))\{ 0 } dsdx 

ill a\Vx(t)\^ + + fo-e(a) “ - 


2 ile,aa I ,2 


ile 


dsdx 


J D 


Lemma 4.3. 0 < —(To,e < 0(e^/^). 


Proof. First, cre(a) is the principal eigenvalue (with principal eigenfunction (p{x, a) = 
iieix,a)) of 

—aA^ — + h^{x)4> = (T0 in fl = Z) X (a,a) 

subject to homogeneous Neumann boundary condition, with a variation charac¬ 
terization analogous to ()4.4p and (j4.5l) . Since the integrand in (14.51) is monotone 
increasing in a < a < a, cre(a) is necessarily less than the principal eigenvalue of 
(lOl) . which is zero. This proves cro,e = cre(a) < 0. 

For the upper estimate, we use a test function 4>{x,s) = ri{s) for (14.91) . where 
T] : [0, oo) —>■ [0,1] satisfies 

r]{s) >0 for 0 < s < 1, and r]{s) = 0 for s > 2. 
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Then upon using S/xV = 0, (I4.2I1 . and 



7?"(0) d 
2 da 



= 0 , 


(since ri{s^) = 0 for e small, and by normalization we have a) dx = 1 for 

all a), we obtain from (SH) that 


-o-o,€ < 


Ip fo V'e + (gi,e£^^^s + 0 ' 2 ,£e^^^s^ + 0 (e^))? 7 ^] dsdx 

Id fo i’Wdsdx 


The conclusion follows from Lemma |4Tjiv) . 


□ 


5. Uniform limit of Uj. 

In this section, we show that converges to 9a in C{D). In particular, hg —>■ 
0a—m\TL C{D). 

Recall that is defined in (14.71) . 

Lemma 5.1. For all /3 > 0, there exists C > 0 independent of e, such that 
We(x, s) < Ce~^e~^^ for all x € D and 0 < s < Se, 
where s^ = {a — ocfje^l'^. 

Proof. First, we derive a rough upper bound of We from Lemma [3.31 
Claim 5.2. There exists C > 0 such that 

sup We < Ce”^. 

D X (a ,a) 


By definition, supwe < (supue)/(inf'0^), and the claim follows from the upper 
bound of Ue (Lemma 13.31) and Lemma I4.11 ivl . 

Next, we construct a supersolution to prove the exponential decay. 


Claim 5.3. For each j3 >Q, there exists sq > 0 independent of e such that 

O'e(a) — cre(a) O’e(o) i/o'4’t,aa(x, o) ^ q 2 f n ^ t , 2/3 —1 

- 273 -UTTs-^ u. ( _^ foT^ all a € [a + € ' so,aJ- 


,2/3 


i>eix,a) 


To see the claim, we note that since cig is monotone increasing in a (Proposition 
IB.lI iiiU. for a = a + and s > sq, 


aeia)-ae{a) cre{a) ^ ae{a +e'^/^so) - cre{a) ao,e 
g2/3 g2/3 — g2/3 ”*”^2/3' 


By (14.21) and Lemma fd.ll il. 

,..rMa + e^/^so)-ae{a) 


lim inf ■ 
e-J-O 


,2/3 


> (lim inf cti (:)so- 
£->■0 


Taking also Lemma [4.31 and Corollary 14.21 into account, we conclude that for a = 
a + and s > so, 


lim inf 
£->■0 


ae{a)-Ue{a) CFeja) 4 / 3 /(’e,aa(a^, a;) 

^ ijje{x,a) 


> (liminf cri,e)so — C. 


£->■0 


Since liminfj-^o > 0 by Lemma [CT hil . Claim lOl holds by choosing sq large. 
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Claim 5.4. For each /3 > 0, there exists sq > 0 independent of e and a supersolution 

( 


W(x, s) := 

y 0 < s < SQ 

defined on D x (so,Se) such that 


sup We I [exp(-/3(s - So)) + exp(^(s - (3/2)se))], 

X ^ D 


(5.1) 


+ W > 0 


dW 

dn 


= 0 


W{x, So) > We{x, So) 

tpt{x,a) 

where Se = (a — . 


^ ll"s(a;,Se) >-e2/3%^^lF(a;,Se) 


in D X (so,Se), 

on dD X (so, Se), 
for X G D, 
for X G D, 


To show the differential inequality, note that the term involving derivatives in x 
vanishes, and that by Claim [57^ and Corollary 14.21 

> + o(l)W, + 2I3^W 

= (-/S^ + o(l)/3 + 2P‘^)W > 0. 

It remains to check the boundary condition on D x {se}, as the rest follows by 
definition. Note that s^ —>■ oo as e —>■ oo, so that exp(—/3(se — so)) exp(—/3se/2). 
We therefore have 

Ws{x,Se) ^ ^ -P exp(-^(se - So)) + 13 exp(-/3se/2) 2 / 3 ^ 

TF(x, Se) ^ ifeix,a) exp(-/3(se - So)) + exp(-/3se/2) 

which converges to a positive constant /3 uniformly for x G D. This proves the 
claim. 

Now, we claim that 

(5.2) Wf^<W inT>x(0,Se). 

By definition, it is easy to see that We < IT in H x [0, so]. That the inequality holds 
in X (so,Se) is due to the fact that IT is a strict positive supersolution of the 
linear problem (EB with homogeneous Dirichlet boundary condition on I? x {so}, 
and Neumann condition on the remaining boundary portions. Standard maximum 
principle applies and shows that the quotient is non-negative. (See, e.g. [Sj 

p. 48].) 

Finally, we obtain Lemma [5.1 1 bv combining Claim [5^ and (15.21) . □ 


Lemma 5.5. Let he given by (|3.2I) . then 

sup |ue(a;) — aue(a:)| —>■ 0. 

x^D 

















AN INTEGRO-PDE MODEL 


13 


Proof. Given e, take S = \/eSi, where (5i is given by Lemma [3.II 


e(x) -aue{x)\ = 


pct 

/ {a — a)Ui{x,a) da 

J Ot 



pct-^-S 


POL 

< d 

/ Uf^{x,a)da 

J 

{a — a) 

/ Ue(x,a)da 
J a+(5 


pCt 

< Sue{x) + C / exp(—/3(a — Q:)/e^^^) da 

J a+i5 

< -v/e + o(l) 


where we have used Lemma I4.1f ivi and Lemma 15.11 in the second last inequality. 
This proves the lemma. □ 


Proposition 5.6. Ue —>■ Oa in C{D), where 6a is the unique positive solution of 
(1^ . In particular, h^^9a — rn in C{D). 


Proof. By Remark 13.21 and Lemma [?751 we deduce that up to a subsequence Cj —>■ 0, 
both Uf_. and Ug/a converges uniformly in D to some uq G W'^''p{D). We claim that 
uq is a (strong and therefore classical) solution of (12.311 . i.e. for each z{x) G C°°{D), 


(5.3) a f IuqAxZ + UQ(m — uo)z] dx — a f uq-;:— dx = 0. 

Jd JdD 9n 

To show ()5.3|) . multiply (13.411 by a test function z{x) and integrate by parts, using 
the Neumann boundary condition of Ve, we obtain 

f [v^A,^z + (m - Jz] dx - f 
J D Jdl 


Ue— dx = 0. 
I dD Dn 


Then one can pass to the limit to obtain (15.311 by invoking Lemma l5.5l By the lower 
estimate in Lemma (3.11 there exists di > 0 such that U[){x) > di for all x G D. 
Hence uq is the unique positive solution of (12.3L i.e. Uej 9a in C{D). Since the 
limit is independent of subsequences, we deduce that that Ue —>■ as e —>■ 0 (not 
just along subsequences Cj —>■ 0). This proves the proposition. □ 


Corollary 5.7. Let a^^a) and il)e{x,a) he the principal eigenvalue and eigenfunc¬ 
tion of and let a*(a) and ip*{x, a) he those of (j2.4ll . Then as e —>■ 0, Ue —>■ a* 

in C'^([a, a]) for all k, and ^e(',a) —>■ '!/>*(•, a) in C^{[a,a\,W'^''P{D)). In particu¬ 
lar, 

ri/T* 

(5.4) (To e ctq := O'* (a) and (Ti ^ crl := ——(a) > 0. 

’ ’ aa 

Proof. Since now h,: ^ ho = 9a — rn in L°°[D), the corollary follows from Proposi¬ 
tion |BTjii). Since ho = 9a — m is non-constant fLemma I3.4|l . Proposition IB . 1 l ivl 
asserts that cr^ > 0. □ 


6. Convergence of 

Let be given by 621), define the normalized version rcg 


by 


We{x,s) 


f ip^idx 
V/nV'>? dsdx J 


We (x, s) on I? X [0, Sj] 


We(x, S) 
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SO that 

( 6 . 1 ) / / ip1{x,a +e^^^s)w1{x,s)dsdx = f 6\dx>{). 

JD Jo JD ~ 

Proposition 6.1. (i) ||t(ie||Lri(Dx(o,Se)) *5 bounded uniformly in e. 

(ii) For any P > Q, there exists Ci > 0 such that for all e sufficiently small, 

We{x, s) < Cie~^‘^ 

for all 0 < s < s,^. 

(iii) e —>■ 0, We{x, s) converges locally uniformly in Dx [0, -l-oo) to the unique 
positive solution of the problem 

fjss + (do - cr^sjf = 0 for s > 0, 

7)s(0) = 0 = 77(-too), f^fj‘^ds = l, 

where do = with ct* given by (lOl) and Aq being the absolute 

value of the first negative root of the derivative of the Airy function. 


( 6 . 2 ) 


Since ? 7 (-|-oo) = 0 and We(x, s) —>■ 0 as s —>■ -|-oo uniformly in x € D, we have in 
fact proved the following. 

Corollary 6.2. ||w<:(a;, s) - i?(s)||L=o(Dx(o.Se)) -t 0 as e -)> 0. 

Proof. By Lemma 14.31 and the fact that r&e is a minimizer of (HU), we obtain 




at -b |(We)sP + 


ID Jo 

+ f [tpe.ai’ew'^] I’Lq dx < C 


cre(a) - cre(a) 


Id 


Claim 6.3. 


£2/3 

V'e Wg dsdx ) = C. 


0(e4/3) 




dsdx 


ID 




D Jo 

< C 


D Jo 


ile (|V. 


We 


+ We) dsdx. 


To prove Claim [631 we apply the Trace Theorem, so that there is C > 0 inde¬ 
pendent of e such that 


ID 


[l/’e.aV'. 


-^e\;=o 


< c 


< c 


D Jo 


D Jo 


(iVajWeP -b \We,a\‘^ + w'^) dsdx 
(|VxWe|^ -b \WeA^ + wf) dsdx, 


where we have used \\'4’e,a'4’e\\L°°{Dx{a,a)) ^ C lCorollarv l5.7p for the first inequality, 
and Lemma Id.ll ivl for the second inequality. 

From Claim 16.31 the normalization m, the estimate in the beginning of the 
proof, and the monotonicity of ae(a) in a (Proposition [BTjiv)), we have 


(l-Ce^/S) 


re 


ae 


-2/3|V, 


We 




I'^+We 


Id Jo 

By Lemma inT ivl. we deduce 


< C 


Ipewl dsdx = C. 


D Jo 


(6.3) 


D Jo 


z -b |(We)sP + 




which implies our assertion (i). Passing to a subsequence, We{x, s) converges weakly 
in HlrD X [0,oo)) to some function fj. of x. Moreover, as \'\/xWe'\^ dsdx < 

Ce^!^, it follows that = 0 a.e.. 
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We outline the rest of the proof of Proposition l6.ll First we will show (iii) except 
the normalization condition 

pOC 

(6.4) / jf ds = 1. 

Jo 

Second, we will show the estimate (ii). Finally we will use (ii) to derive ()6.4I1 from 
(EH), which completes the proof of (iii). 

We claim that fj must satisfy the equation in EH- To see this claim, note that 
the equation for die is 


(6.5) 


0 = -- 






We argue via the weak formulation. 


Claim 6.4. There exists a constant uq such that for each test function z(s) that is 
compactly supported in [0, oo), 

[-Zsfjs + (do - cr*s)z77] ds = 0. 

In particular, fj satisfies the equation fjss + (do — cr*s)Tj = 0 on (0, oo) and fjs(O) = 0. 



Multiplying (16.51) by a test function z = z(s), and integrating over x € D, we 
see that the term involving derivatives in x vanishes (by the Neumann boundary 
condition = 0), and obtain 
( 6 . 6 ) 

'Ip'^We dx. 


0 = -z [V’e (r5e)s]s dx + z / 

J D J D 


ae(a) - aeia) (Te(a) AI3i’e,a. 
g2/3 + g2/3 ^ 

Next, integrate the first term of ()6.6I) over s G [0, s^], we see that 

- Z [lj}‘^(We)s]sdxds 


Id 


Zsl(j'^{We)s dxds — 


lo JD 


Z / ’4fl(We)s dx 


ID 


J s=0 


Zs%lfl(we)s dxds + 


/O JD 


Z / 4’e'4’€,ayje dx 


ID 


J s=0 


where we have used the boundary condition (w^)s = on Dx {0, Sg}. 

Since z(s) has compact support in [0,oo), the boundary term evaluated at s = s,; 
vanishes, and the remaining boundary term is of order 0(e^/^) (since We is bounded 
in H^(D X (0, Se)) by assertion (i), and hence bounded in Lf(D x {0}) by the Trace 
Theorem). Hence, we have 


(6.7) 


z / [il}‘^(w,,)s]sdxds = 


D 


lo JD 


Zs'ipl(We)s dxds + o(l). 


Also, in the support of z(s), ('0£)^(x,a + e^/^s) —>■ (V’*)^(x, a) uniformly, so we may 
use (16.7p to integrate (16.611 over s S [0, Se] and pass to the limit to get 


( 6 . 8 ) 


0 = 


('ip* Y(x,a) dx 


’D 


poo poo 

/ Zsfjsds-\- / {a^s — ao)zf]ds 

Jo Jo 
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where we have used Corollarv l5.7l and that ao is a subsequential limit of —CTe(a)/e^/^ 
(see also Lemma lOl) . This proves Claim 15^ Next, we claim that 


(6.9) 


fj^ ds < + 00 . 


Jo 

Notice that by normalization of (see (16.11) 1. and the uniform (in e) positive up¬ 
per/lower bound of V'e (Lemma Id.lf ivll. there exists a fixed constant Co such 
that for each M > 0, and for all e > 0 sufficiently small, w1 dxds < 

JoW^dxds < Co- Letting e —>• 0, if ds < Cq for all M > 0. i.e. (16.91) 
holds. 


Claim 6.5. fj is a positive solution that satisfies (16.2p with condition (16.41) being 
replaced by (ESI). 


By Claim 113 ij satisfies 

(6.10) rjss + (ao - Cis)?) = 0, ?) > 0 on [0, +oo), and 7)s(0) = 0. 

It remains to show that p^+oo) — 0, and that the subsequential limit Oq must 
be determined by do of the proposition. By (16.101) and i] > 0, ijss > 0 for all s 
sufficiently large. Hence 7)(+oo) exists in [0,+oo]. By (16.91) . ri{+oo) = 0. 

Hence, ?) is a constant multiple of Airy((cr*)^/^s — Aq), where Ao is the absolute 
value of the first negative root of the derivative of the Airy function Airy(a;). In 
particular the subsequential limit do given in (|6.8I) is uniquely determined by do = 
{a*f^^Ao (i.e. the full limit limg_,,o (Te(a)/e^/^ exists). This shows Claim [6+1 To 
finish the proof of (iii) except for (16.41) , it remains to establish the following. 

Claim 6.6. Weix, s) —>■ r)(s) locally uniformly in Dx [0, oo). In particular, for each 
M > 0, ||u)e||L=(_Dx[o,M]) bounded uniformly in e. 


It is enough to show that for each M > 0, 

( 6 . 11 ) 


S^VxeDWeix,s) 

sup ^->• 1 


sG[0,M] ide(x, s) 


For, assuming (16.111) . one can write 


as e \ 0. 


(6.12) i(}e(a;, s) = u;£(a;o, s)(l + (5e(a;, s)) for some ccq G fd, 

where ^^(a:, s) —?> 0 in x [0,+oo)). Now, if we integrate (16.121) over x £ D, 

then 

iLe(s) I We{x,s)dx = wfixQ,s){l + Sfis)), 

where Sfis) —^ 0 in L^^([0,oo)). Since We is bounded in H^{D x (0, s^)), one 
can easily deduce that We{s) G iL;L((0,+oo)) C {[0,+oo)). Therefore, by 
Arzela-Ascoli Theorem, Wfis) and hence wfixo, s) converges to r]{s) in C'ioc([0, oo)). 
Finally, (16.121) implies that wfix,s) —>■ 7)(s) locally uniformly in D x [0,+oo). 

It remains to show (16.lip in a similar fashion as in Claim [23 Assume to the 
contrary that there exists some constant cq > 1, e = —?> 0, —>■ sq < +oo, such 

that 


sup We{x,Sk) > Co inf We{x,Sk)- 

xeD 


(6.13) 
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Similarly as before, we extend We by reflection on Z3 x {0} so that We is defined on 
D X (—Se,Se), and define 


Wk{x,T) 


Wejx^Sk +e'^^^Tl^) 
Snp^'eDWeix',Sk) 


in D X 




{Se 



Recall that Se is defined in (I4.6I1 . By the equation (16.5|) satisfied by We, Wk satisfies 


(6.14) 


in Z) X I — 


-|V, • - {lPlWk,r) 

+i 




,2 be 


^IWk = 0, 


X (^-(s, + Sfc)^,(se-Sfc)^), where 


a = air) = a + 


Sk + e- 


and the boundary conditions 

f 0 X (-(Se + Sfe)^,(Se - Sfc)^) 

[ Wk,T = -^^Wk onZ)x |-(se + Sfc)^,(Se-Sfc)^|. 

Since Sg —>■ +oo as e —>■ 0 and that Sk remains bounded, we see in particular that 
the domain of Wk tends to I? x K. as fc —>■ oo. 


Claim 6.7. For each M > Q, 


a + +e —— 


—>■ 0 as e —>■ 0, uniformly for 


T G [-M,M]. 


By Lemma m, CTj are bounded in C^{[a,a]) uniformly in e. Hence we may 
write 


ae ( a + SkCS + e-^ ) 

< \creia)\ + C 

e^l^SkFe^ 

V v^y 




and conclude that + SfeC^ + goes to zero by Lemmaand bounded¬ 
ness of Sfc, T. This proves Claim 1^771 

Since the coefficients of (16.141) are bounded in x R) uniformly in k, Har- 

nack inequality, and the normalization condition sup^.^^) Hfc(a:,0) = 1 implies that 
Wk are bounded in x R) uniformly in k. Hence we may apply elliptic LP es¬ 

timates similarly as in Claim [3^ to conclude that a subsequence of Wk converges in 
L^^{D X R) to a positive solution of ('0o(a;, a))“^Va; • {'f’Q{x,(^VxW) + Wrr — 0 in 
H X R. (Here we used Claim IHTTl l Now, we apply the following Liouville Theorem, 
whose proof is exactly analogous to Proposition lC.il and is skipped. 


Proposition 6.8. Let f^{x) be a smooth positive function defined in D, then every 
positive solution W to x • {'4’^'^xW) + Wtt =0 in D x R, subject to Neumann 

boundary condition on dD x R, is necessarily a constant. 


So that by normalization sup,j.g£, H4(a;, 0) = 1, H4(x,0) —>■ 1 uniformly in D. 
This contradicts (j6.13l) and proves ()6.1ip . This establishes Claim [RH Claims l6^ 
and 16.51 establish assertion (iii) except for condition ()6.4I1 . 

Next, we proceed to show the estimate in (ii). By the preceding argument in the 
proof of Lemma [531 specifically the construction of supersolution W in Claim [531 
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sup I [exp(-/3(s - So)) + exp(/3(s - (3/2)se))] 

X ^ D 


we can show that for all /3 > 0, there exists so > 0 such that 

/ 

w{x, s) < 

y 0 < s < SQ 

for X & D and s G [so, Se). Then (ii) follows from Claim [6761 as the expression inside 
paranthesis is bounded uniformly in e. We do not repeat the details. 

For (iii), it remains to show (16.41) . By assertion (ii), and that 

(6.15) tpeix, a + ^'tjj*{x, a) and Weix, s) ^ fi{s) 

in X [0, oo)) (by Lemma HTTT ivi and Claim I57BI resn.i. we may pass to the 

limit in (16.11) to obtain 


ID 


6^dx= / tp'^{x,a + e'^^^s)w'^{x,s)dsdx-^ / {'ip*Y{x,^dx / if ds 


D Jo 


ID 


Upon noting that (see Definition USTii)) 

(6.16) 'll)* {x, a) = Ogfx) in D, 

the proof is completed. 

7. Proof of Theorem 12.31 

Proof of r/teorem 12.31 First, we note that by Proposition 15.61 
(7.1) 


□ 


ipeWtdsdx = / / ufx,a)dadx= / u^dx ^ / 9a dx 

D«^0 J D J ^ J D J D 

as e —>■ 0. Furthermore, by (I6.15|) . (16.161) and the estimate of Proposition 16.If iil. 


(7.2) 


D Jo 


ipeWe dsdx—>■ / il)*{x,a)dx / fj{s) ds = / 9gAx)dx / fj{s)ds. 
J D Jo JD JO 


By the definition of We and there is a function r(e) such that 

(7.3) r(;e(a;, s) = F(e)r(;e(a;, s). 

By jlU and (177^ . we have 

(7.4) lime2/^r(e) = ?7ds^ . 

Hence, by (17.31) and Corollary 16.21 

(f/^'wfx, {a - a)/e2/3) “ ^ 9 


0 . 


L”(n) 


By the fact that fjds) ^ ii{s) = ri*{s) where rj* is given in Definition iii), 
we also have 


t^^^ivfx, (a — a)/e^^^) ~ 9* 
Using Lemma l4.1f ivl. we have 


a — a 

f2/3 


0. 


(7.5) 




fx, a) - ipeix, a)r]* 


L“(n) 


L“(n) 


0 . 
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By the fact that r]*{s) < Ce for some C,/3 > 0, (16.151) and (j6.16|) . we have 


(7.6) 




0 . 


L“(n) 


And p.7p follows from (17.51) and (j7.6l) . 


□ 
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Appendix A. Existence Results 

In this section we show the existence of positive solution to (EH). For this 
purpose, we fix positive parameters e and a> a, and denote (in this section only) 
the principal eigenvalue and eigenfunction of the following problem by /ii and (j)i. 

( aA^cj) + e^(j)aa + m{x)(p + iKp = 0 in := D x (a, a), 

§^=0 ondDx{a,a), 

4>a = 0 on D X {a, a}. 

Theorem A.l. ///ii > 0, then the equation EH has no positive steady-state. If 
pLi < 0, then the equation EH has at least one positive steady-state. 


Proof. First, we prove the non-existence result. Suppose > 0 and let it be a 
non-negative solution of EH- Multiply EH by the principal eigenfunction (fi of 
and integrate by parts, we obtain 


Ml 


dadx -\- 


uu(pi dadx = 0. 


Since /ii > 0, both terms are non-negative, and both must be identically zero. i.e. 

It = 0. 

For the existence result, we consider, for r G [0,1], the positive solutions of 

. . J aAu -\- e^{u)aa + {m{x) — tu — {1 — t)u)u = 0 in £> x (a, d), 

\ ^ = 0 on dD X (a,d), {u)a =0 in 71 x {a, a}. 


Here we recall that u = u(x, a) da. It remains to show the following claim, 
from which existence of positive solution to EH follows by a standard topological 
degree argument, as the existence of a unique, linearly stable positive solution to 
(IA.2|) when r = 0 is standard. 


Claim A.2. For some (5 > 0 independent of t G [0,1], any positive solution u of 
(Q satisfies 


^ < IhllLi(n) < l/<5- 
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For the upper bound, one can integrate (IA.2I1 over V, to get 


/ umdx = T / dx + {1 — t) 
J D J D 


dadx 


D J OL 


> T 


1 — T 


dx 


JD 


> Co / dx 
JD 


> 


Co 

\D\ 


' D 


udx\ = 


Li(n)5 


from which the upper bound follows. 

For the lower bound, let u = v4>i^ where > 0 is the principal eigenfunction 
corresponding to the principal eigenvalue < 0 of (lA.ip . Moreover, if we normalize 
(j)i by /q (/)f = 1, then sup^ (/)i and infn are fixed positive constants independent 
of r, as (jA.ip is independent of t. The equation for v can be written as 

X ■ (<^fVa;u) + e^{(t)\va)a + “ Tu^ “ (1 “ t)u) = 0 in O, 


dadx > 0. 


= 0 on dD X (a, a), Ua = 0 on Z) x {a, a}. 
Hence, if we divide by v and integrate by parts, we have 

2Q!|Va;Up + 


/ (j)'f{^i+TU+{l — T)u)dadx= / 4>l- 

Jci Jn 


Hence we have 


^sup(/)i^ [T(a - a) + (1 - r)]||u||Li(o) > -//ly (jif dadx = -fii > 0. 


Since /ii and supf; (j)i are independent of r, we have 

-Ml 


klUqn) > 


(supn (^i)2 [T{a - a) + (1 - r)]' 


Since the latter term is bounded from below uniformly in r G [0,1], the claim is 
proved. □ 

Corollary A.3. If m{x) dx > 0, then for any e > 0, (12.11) has at least one 
positive solution. 

Proof. Divide the equation (lA.ll) by the principal eigenfunction (j)i and integrate 
by parts over fl, we get 


a\Vx4>i\'^ + 


dadx + / / ( 771 ( 0 ;) +/Ti) dadx = 0. 

J D J a 


Hence for all e > 0, 


Ml < “71=^ / 'm{x)dx < 0, 

1^1 JD 

and the existence of positive solution of (12.11) follows from Theorem lA.il 


□ 
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Appendix B. Eigenvalue problems with diffusion parameter a and 

WEIGHT FUNCTION h{x) 

For each a > 0 and h G L°°{D), let Ai = Ai(a, h) G R and ^p{x) = /i) 

be the normalized principal eigenvalue and principal eigenfunction of the following 
problem. 

J -aA^ip + h(p = Xip in D, 

11^ = 0 ondD, f^ip^dx = l. 

We shall state and prove a number of properties of Ai and (pi, and its dependence 
on the parameters a and /i, some of which is folklore among specialists. 


Proposition B.l. (i) For each a > 0 and h G L°°{D), the problem (IB. II) has 
a principal eigenvalue Ai which is simple, and the corresponding eigenfunc¬ 
tion ipi can be chosen positive and uniguely determined by the constraint 
foFldx = 1 . 

(ii) For each p > 1, the mapping {a,h) i—>■ (Ai,(^i(-)) is smooth from R+ x 

L°°{D) toRxW^P{D), where W^P{D) = {(!)£ M = 0 on dD}. 

(iii) If h G L°°{D) is non-constant, then ^^(a,h) > 0 for all a > 0. 


Proof. Part (i) is well-known. See, e.g. [141 Section 8.12]. In particular, the 
principal eigenvalue is given by the variational characterization 


(B.2) 


, , ■ r + h(p^) dx 

Xi{a,h)= inf -p-—-. 

vpgC1(d)\{0} JdF dx 


Fix p > N {N being the dimension of D). Consider the following mapping F : 
X R X R+ X L°°{D) -g LP(L>) x R, given by 

F{ip, A, a, h) = {aA,cP - hp-\- Xp, / dx - 1), 

J D 

Then for each a > 0 and h G L°°{D), the principal eigenpair a, h), Ai(a, h)) 

of (jB.lIl satisfies 


(B.3) F{pi{-;a,h),Xi{a,h),a,h) = (0,0). 


Assertion (ii) follows from the following claim, in view of the Implicit Function 
Theorem and the smooth dependence of the operator F on a and h. 


Claim B.2. For each fixed a > 0 and h G L°°{D), 

D^,p^x)F{pi{--,a,h),Xi{a,h),a,h) : IF^^ x R L^{D)xR 
is a bijection. 

We shall follow the arguments in the proof of nni Lemma 2.1]. Suppose for some 
($,t) G X R, D(^,^^x)F{pi{--,a,h),Xi{a,h),a,h)[^,i\ = (0,0), i.e. 

5 $ 

(B.4) aAa;$ —-b Ai$ + = 0 in £>, and -— = 0 on dD, 

an 

and 

(B.5) 2 / ^pxdx = Q, 

JD 

where Ai = X{a,h) and pi = p{x;a,h). The Fredholm alternative implies that 
Jjjtpfdx = 0, i.e. t = 0. Hence $ = cpi for some constant c (as Ai = Ai(a, h) 
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is a simple eigenvalue). But then (IB.5I) implies that c = 0. This shows that the 
kernel of •; a, h), Ai(a,/i), a,/i) is trivial. Now let {f,q) G LP{D) xR 

be given, we need to solve for ($,t) in 

5 $ 

(B.6) aAa;$ —+ Ai<i>-b = / in £>, and = 0 on dD, 

an 

and 


(B.7) 2 / dx = q. 

JD 

Set t = (/^ fifi dx)/{Jjj dx), then Jjj{f — dx = 0, so (|B.6(1 has solution of 

the form where G W^^{D) is unique and satisfies dx = 0. 

Finally, if we set s = q/{2 ip\ dx) then {^,t) solves (IB.61) and (IB.7|) . This proves 
Claim IB?^ which implies assertion (ii). 

For (iii), we differentiate (IB.Ill with respect to a, 


r -aA^^ + h^-X,^ = in D 

I =0 on dD, and dx = 0 


Multiply (IB. 81) by (pi and integrate by parts, we have fjj (pf dx = \Vdx. 
Since h{x) is non-constant in x, = ipi{-; a, h) is non-constant in x and this 
implies that > 0. This proves (iii). □ 


First, we show that Ai and gji are continuous with respect to the weak topology 
of M+ X np>iLP{D). 

Lemma B.3. Let Ai(q;, h) and ipi{- ;a,h) be the principal eigenpair of (IB.1|1 . 

(i) For each p > 1, there exists Cq = Cq{p, M,a,a,dD) such that 

\Xi{a,h)\ + \\pi{- ■,a,h)\\^^2.p^D) < C'g 

provided a G \a,a] and < M. 

(ii) If aj —>■ oq G [a, a], supj>Q ||/ij 11^00(23) < +oo and hj ho in Lp{D) for 
all p > 1, then as j —>■ oo, Xi{aj,hj) —>■ Ai(q;o,/io) ipi{ - ;aj,hj) 

■ ;ao; ho) weakly in W^'P{D) for all p > 1. 


Proof. By (IB.21) . Ai := Ai(q;, h) forms a bounded sequence in [—||/i|lioo(£)), 

The LP estimate (for p > N) applied to (IB.II) and interpolation inequality together 
imply ^ 

\Wl\\w^<P(D) < C\\ipi\\LP{D) < 2\\Fi\\w^^p{D) + C'||i^i||l2(D), 

where C is a generic constant, depending on \\h\\L^(^jopa,a and the domain D that 
changes from line to line. This proves (i). 

For (ii), let aj —>■ oq G [a, a] and hj be a uniformly bounded sequence in (D) 
and hj ho weakly in LP{D). Denote Aij- = Xi{aj,hj) and = (pi{ - ;aj,hj). 
By assertion (i), there are subsequences Aij/ and pij' such that Xi{aj,hj) —>■ A 
and pi{- ] aj,hj) ^ (fi weakly in W'^’P{D), for some A G R and p G W^’P{D). Take 
a = aj', h = hj! in (IB.II) . and pass to the weak limit f ^ oo, we deduce 


—aoAxP + hop = Xp in D, 
1^=0 ondD, fjy<p^dx=l. 


Hence (p, X) is an eigenpair of (IB. Ill when a = ao, h = ho and such that (p > 0. 
Moreover, <p is non-trivial, as Jjjip^dx = 1. By uniqueness of principal eigenpair. 
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it follows that A = Ai(ao,ho) and (f = (pi{- ;ao,ho). Since the limit is indepen¬ 
dent of subsequence, we deduce that the full sequence Xi{aj,hj) —)> X{ao,ho) and 
ipi{ -; aj,hj) ^ (/?i( •; Oq, ^o) weakly in W^’P{D). This proves the assertion (ii). □ 


Next, we show the following uniform estimate of 
Lemma B.4. There exists C 2 = C 2 {M,a,a, D) such that for any a € [a, a] and 

II^IIl“(d) < if 

D(,p,\)F{(pi{-;a,h),Xi{a,h),a,h)[<^,t] = {f,q), 
i.e. dEH) and (EJl) hold with Ai = A(a, h) and Lpi = Lpi{x] a, h), then 
(B.9) |t| -I- ||$||vv2.p(D) < C2{\q\ + ||/||lp(d))- 


Proof. Let M > 0 be given. Suppose to the contrary that there are aj € [a, a ], hj , 
$j, tj , qj , fj such that 


(B.IO) sup < M, \tj\ F ||$j||w2,p(d) -t cc, \qj\ + \\fj\\LP{D) < 1- 

3 

Without loss of generality, we may assume aj —>■ oq S [a, a] and for some ho € 
L°°{D), hj ho weakly in U'{D). Denote 

Ai,j = Ai(Q;j, hj), and ipij = ipi{-]aj,hj) for j G N U {0}. 

The above arguments ensure that 

^j = ^f+qj/{2 ipljdx)ipij, and tj = / fjipijdx/ Lp\ j dx 

JD ’ JD JD 

where is the unique solution of (IB.61) subject to the constraint dx = 0. 

By the normalization (pf ^ dx = 1., we have 


(B.ll) 

^ JD 

Since we have shown that |Aij| and remain bounded uniformly in 

j, (ITOl) and (IB.Ill) imply that \\^j-\\^r 2 ,p(^Dj —>■ 00 . Apply estimate to the 
equation of which is 


(B.12) 


- hj<^f + Aij$-L = fj - (D dx)(pij in D, 

= 0 on dD, and dx = 0. 


Using the boundedness of (pij in W^’^{D) and hence in L°°{D), we have 


\^f\\w^,p{D) < C 


^j‘llLP(D) + ll/j - fjPl.j dx^ PljWLPiD) 


< C'dl^j IU“(D) + |l/j||LP(D))- 

Hence we must have ||$j‘||L°p(D) —>■ 00 as j —>■ 00 . Define 4>j := 
then $j satisfies 


f - hj^j + Aij$j = fj in D, 

1 Id dx = 0, and sup^, l>j = 1, 

where fj = [fj - {Jjj fjTi.j dx)(pij]/\\^f\\L^(^D) converges to zero in LP{D) as 
j —>■ 00 . By PP estimates, $j is bounded uniformly in W^’P{D). Hence, there is 
a subsequence $j/ that converges, weakly in W^’P{D) and strongly in C^{D), to 
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some function $o- By normalization sup^) <i>o = lim^v (^sup^, j = 1. Moreover, 
$0 satisfies (using Lemma |R3Kii)) 

f ao^x^o — ^ 0^0 + Ai(ao: ^o)®o = 0 in £>, 

1^=0 Id ^oV^lo dx = 0. 

Since $o is non-negative, Proposition IB.ir il implies <f>o = c(pi{- ;ao,ho) = C(pifi 
but the integral constraint implies that c = 0. i.e. <i>o = 0- This is a contradiction 
to sup^, l>o = 1. This proves (jB.9|) . □ 


Proposition B.5. Let Ai(a, h) and ipi{- ]a,h) be the prineipal eigenpair of (|B.1|) . 
(i) For each k, there exists C'f. = C'^.(M,a,a.,D) such that 


(B.13) 


k 

E 

4=0 





+ 




<c; 


IV2,p(£)) 


provided a € [a, a] and ||h||ioo(-£,) < M. 

(ii) //supj>Q ||hj < - 1-00 and hj ho in LP{D) for all p > 1, then for 

each k > 0, 

Qk Qk 

^^Ai(-,hj)-)> ^^Ai(-,ho) in Cioc{[0,oo)) as j ^ oo. 

Moreover, given k > {), p > 1, and sequence aj —>■ ao > 0; 

■aj,hj) -- ■,ao,ho) weakly in W‘^'P{D) as j ^ oo. 


Proof. Assertions (i) and (ii) for the case k = 0 are exactly Lemma FB.31 We first 
prove assertion (i) for fc = 1, differentiate the relation (IB.3|) with respect to a, 


(B.14) 




L^,(-,c,.h),Lxoc,,h) 


= -D„F. 


where the partial derivatives of F are evaluated at (}pi{ ■ ;a,h), Ai(a, h), a, h). By 
M . we may write 

((^;,a;) = (z?(^.;,)F)-i[-D„i^] = (i?(^,;,)F)-i(-A,<^i,o) 
and deduce by Lemma FB.41 that 


l-^ll + II‘PiI|iU2.p(D) < C2\\^x‘Pi\\lv{D) < C2\\'Pi\\w^’P(D) < C. 

i.e. assertion (i) holds for fc = 1. We argue inductively ior k > 1. Suppose (i) holds 
for k = K — 1. We can write 


(B.15) 

where 




- gK qK 

da^ ^i(«i h) 


= ^k{q:, h) 


(B.16) SK{a,h):= 


gK-i ^-1 

-K ^ - Yi 

k=l 


da^ 


K 

k 


gK 


-k 


gk 

da^ ^da^~^ 


Pi 


By the form of we can deduce the following result. 


Claim B.6. ||5K(a,/ i)||l~(d) < 



doL^ 


pi 


IU2.p(D) 
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By the induction assumption (i.e. (i) holds ioi k = K—1) we have ||lp(_d) < 
C{M,a,a,D). Hence we may apply Claim [BH] to (IB.15I) to conclude the assertion 
(i) for the case K. This induction argument proves (i). 

By Lemma |B.3f iil. it remains to prove assertion (ii) for case fc > 1. Let aj —>■ 
ao G [«, a] and hj be a uniformly bounded sequence in L°°{D) and hj ho weakly 
in LP{D). Denote Aij- = Xi{aj, hj) and (pij = ipi{ -; aj, hj). By assertion (i), there 
are subsequences Aij' and tpij/ such that for all k > 0, -^j^Xi(aj> ,hji) —>■ Xk and 
- ;aji,hji) (pk weakly in W‘^’P{D), for some A^ S K and p>k € W'^'^{D). 
Passing to the limit in (IB.14I) . we deduce that 


(B.17) 




<f>i,Xi 


= -DryF. 


where the partial derivatives of F are evaluated at (p>i{ ■; ap, ho), Ai(aoi ho), ap, ho). 
Since we also have 


(B.18) 


D{v,x)P 


—(^i( •; ao, ho), —Ai(ao, ho) 


= -DaF, 


where the partial derivatives of F are evaluated at {p>i{ ■; ao, ho), Ai(ao, ho), ao, ho), 
we may invert in both (IB.171) and (IB.18I1 . and conclude that 


‘Pi — ■ ’ “O; ho) 


d 

and Ai = —Ai(ao,/io). 


Since the limit is determined independent of the subsequence, we conclude assertion 
(ii) for the case k = 1. 

Again, we may argue inductively for fc > 1. Suppose (ii) is proved for k = 
1,..., K — 1. The following can be easily observed from (|B.16p . 


Claim B.7. If assertion (ii) holds for k = 1,..., K — 1, then 

^K{oLj,hj) ^k{olo, ho) 
weakly in Lp{D), where is defined in (jB.16l) . 

Based on Claim I bTtI and the assertion (ii) for the cases fc = l,...,Ar—1, we may 
pass to the limit in (IB.lsp . Together with the uniform boundedness of '■ 

LP{D) —>• W'^'P{D) (Lemma IB.41) . this implies ■^jprXi{aj,hj) —>■ Ai(ap,/ip) and 

dK 

-^^ipi{--,aj,hj) p)i{--,ao,ho) in W‘^’^{D). 

Thus assertion (ii) follows by induction on k. □ 


Appendix C. Liouville Theorem for Positive Harmonic Functions in 

Cylinder Domain 

We give a proof of the Liouville-type theorem for positive harmonic functions in 
cylinder domains, since we cannot locate a proper reference for this result. 

Proposition C.l. Let k G N, D be a bounded smooth domain in M.^ and u be 
a non-negative harmonic function on LI := D x C so that ^ = 0 on 

dD X Then u is necessarily a constant. 
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Proof. Let x G D, y G and let u{x, y) be a non-negative harmonic function on 
n = £) X subject to Neumann boundary condition on dD x By subtracting 
a positive constant from u, we may assume that info u = 0. 

Harnack inequality says that there is a constant C > 1 such that for all y' gM.^ , 
we have 

sup u < C inf u. 

a;g_D,|y-y'|<2 x&D,\y-y'\<2 

Define v{y) = jklo u{x',y)dy, then v is a harmonic function on K.* and must 
be equal to a non-negative constant wq- Hence for each y' G M*, there exists x' G D 
such that u{x',y') = vq. It follows that for each y' G 

vo<C inf u{x,y). 

x€D,\y-y'\<2 

Taking infimum in y' G it follows that from info u = 0 that vg = 0. Hence, 

y u(x,y)dx = v(y) =vo = 0 

for all y G i.e. u = 0 in D. □ 
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